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, Using Coherent-State (CS) techniques, we prove a sampling theorem for holomorphic 

^SJ ' functions on the hyperboloid (or its stereographic projection onto the open unit disk Di), 

seen as a homogeneous space of the pseudo-unitary group SU (1,1). We provide a reconstruc- 
(— I ' tion formula for bandlimited functions, through a sinc-type kernel, and a discrete Fourier 

^ |. transform from N samples properly chosen. We also study the case of undersampling of 

"quasi-bandlimited" functions and the conditions under which a partial reconstruction from 
N samples is still possible and the accuracy of the approximation. 
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Abstract 



Introduction 



^ ■ In a previous article [T] , we proved sampling theorems and provided discrete Fourier transforms 

! for holomorphic functions on the Riemann sphere, using the machinery of Spin CS related to 

the special unitary group SU{2), which is the double cover of the group 5*0(3) of motions of 

■ the sphere Here we study similar discretization problems for its noncompact counterpart 
. SO{2, 1) (the group of motions of the Lobachevsky plane) or, more precisely, for its double 

CD I cover SU{1, 1). Both, SU{2) and SU{1, 1), appear as the underlying symmetry groups of many 

physical systems for which they constitute a powerful computational and classification tool. 
In fact. Angular Momentum Theory proves to be essential when studying systems exhibiting 
rotational invariance (isotropy). In the same manner, the representation theory of SU {1^1) 
^ ' or SL(2, M) is useful when dealing with systems bearing conformal invariance, specially in two 

■ dimensions, where this finite-dimensional symmetry can be promoted to an infinite-dimensional 
one (the Virasoro group). In particular, the group 5L(2,M) was used in p] to define wavelets on 
the circle and the real line in a unified way. Furthermore, SU{2) and SU{1, 1) CS, generalizing 
canonical CS of the Heisenberg-Weyl group (Gabor frames), find a great variety of applications, 
mainly in the study of quantum mechanical systems and their classical limit (see e.g. [31 HI [5]). 
For example, ground states of superconductors and superfluids (like Bose-Einstein condensates) 
are coherent states. 

Also, standard Continuous Wavelet Theory (see e.g. [6]) can be seen as a chapter of CS on 
the group of affine transformations (translations and dilations). Here the discretization process 
turns out to be essential for computational applications in, for example, signal processing. These 
results revived interest in the question of discretization and we hope that the establishment of 
new sampling theorems for harmonic analysis on non-Abelian groups and their homogeneous 
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spaces will be of importance for numerical study and simulation of those physical systems bear- 
ing that symmetries. Actually, there are some important general results about sampling and 
efficient computation of Fourier transforms for compact groups (see e.g. [3 [5]). However, a 
comprehensive study of the non-compact case is far more complicated, although some results in 
this direction have been obtained for specific groups (see e.g. [S] for a survey). For instance, 
we would like to point out Ref. [10] for the motion group and its engineering applications [TT] 
(namely in robotics [12]) and p[3| for discrete frames of the Poincare group and its potential 
applications to Relativity Theory. 

This article intends to be a further step in this direction. Completeness criteria for CS 
subsystems related to discrete subgroups of SU{1,1) have been proved using the theory of 
Automorphic Forms (see e.g. [3]). Here we shall follow a different approach. Working in the 
open unit disk Di = SU{1, 1)/C/(1) (as an appropriate realization of the Lobachevsky plane or 
the hyperboloid), we shall choose as sampling points for analytic functions inside Di (carrying a 
unitary irreducible representation of SU (1,1) of Bargmann index s) a set of equally distributed 
points on a circumference of radius r < 1. For bandlimited holomorphic functions on Bi of 
bandlimit M < N and index s, the resolution operator A is diagonal, providing a reconstruction 
formula by means of a (left) pseudoinverse. The Fourier coefficient can be obtained by means 
of the (rescaled) Fourier transform of the data, allowing for a straightforward fast extension 
of the reconstruction algorithm. The reconstruction of arbitrary (band-unlimited) functions is 
not exact for a finite number A^ of samples. However, for fast-decaying, or "quasi-bandlimited" , 
functions it is still possible to give partial reconstruction formulas and to analyze the accuracy of 
the approximation in terms of N, the radius r and the index s, this time through the sampled CS 
overlap (or reproducing kernel) B (see later on Sec. [2]for definitions), which exhibits a "circulant" 
structure and can be easily inverted using the properties of the Rectangular Fourier Matrices 
(RFM) and the theory of Circulant Matrices |14j . This helps us to provide a reconstruction 
formula accomplished through an eigen-decomposition B = !FD!F^^ of B, where J- turns out to 
be the standard discrete Fourier transform matrix. 

The plan of the article is as follows. In order to keep the article as self-contained as possible, 
we shall introduce in the next section general definitions and results about CS and frames based 
on a group G. The standard construction of CS related to the discrete series of representations 
of G = SU{1, 1) is briefly sketched in Sec. El We refer the reader to Refs. P H H [15] for 
more information. In Section U] we provide sampling theorems, discrete Fourier transforms and 
reconstruction formulas for bandlimited holomorphic functions on Di of bandlimit M and index 
s, and discuss the effect of over- and under-sampling. For quasi-bandlimited functions these 
reconstruction formulas are not exact and we analyze the error committed in terms of A, r and 
s. Finally, Sec. [5] is devoted to conclusions and outlook. 

2 Coherent States, Frames and Discretization 

Let us consider a unitary representation [/ of a Lie group G on a Hilbert space (W, (•!•)). 
Consider also the space L?'{G,dg) of square-integrable complex functions ^ on G, where dg = 
d{g'g), yg' G G, stands for the left-invariant Haar measure, which defines the scalar product 




(1) 
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A non-zero function 7 G is called admissible (or a fiducial vector) if T{g) = {U{g)'y\j) £ 
LP'{G, dg), that is, if 



cy= / T{g)T{g)dg= / \{U {g)-im' dg < ^ . (2) 
JG JG 

Let us assume that the representation U is irreducible, and that there exists a function 7 
admissible, then a system of coherent states (CS) of 7i associated to (or indexed by) G is defined 
as the set of functions in the orbit of 7 under G 

lg = U{gh, g£G. (3) 

We can also restrict ourselves to a suitable homogeneous space Q = G/H, for some closed 
subgroup H. Then, the non-zero function 7 is said to be admissible m.od{H, a) (with a : Q ^ G 
a Borel section), and the representation U square integrable mod{H,a), if the condition 

\{U{a{qMij)\^dq<^, G 7^ (4) 

holds, where dq is a measure on Q "projected" from the left-invariant measure dg on the whole 
G. The coherent states indexed by Q are defined as Ja{q) = U{a{q))'y,q G Q, and they form an 
overcomplete set in Ti. 

The condition (j4]) could also be written as an "expectation value" 



< 



/ \{U{a{qM^l^)\^dq={i;\A^\i;) G (5) 

JQ 



where = Jq \ ja{q)) {la{q)\dq is a positive, bounded, invertible operatorQ 

If the operator is also bounded, then the set So- = {|7o-(g)))9 ^ Q} is called a frame 
(see [16] for details on frames), and a tight frame if A^^ is a positive multiple of the identity, 
A^ = XI,\> 0. 

To avoid domain problems in the following, let us assume that 7 generates a frame (i.e., that 
A~^ is bounded). The CS map is defined as the linear map 



T^: n — > L\Q,dq) 



(6) 



Its range L'^{Q, dq) = T^iTL) is complete with respect to the scalar product ($1^)^ = (<I>|T^j4^^T^ ^^)q 
and Ty is unitary from Ti onto L'^{Q,dq). Thus, the inverse map T^^ yields the reconstruction 
formula 



T-^^y= [ ^^(g)A;Sa(g)d'7, ^yeL^y{Q,dq), (7) 



which expands the signal in terms of CS A'^^^^f^^-^ with wavelet coefficients ^'7(0') = [T^tpjiq). 
These formulas acquire a simpler form when A^^ is a multiple of the identity, as is for the case 
considered in this article. 



*In this paper we shall extensively use the Dirac notation in terms of "bra" and "kets" (see e.g. [2] [9]). The 
Dirac notation is justified by the Riesz Representation Theorem, and is valid in more general settings than Hilbert 
spaces of square integrable functions . 
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When it comes to numerical calculations, the integral A^^ = /g |7(T{g))(7(7(q)l'^^ to be 
discretized, which means to restrict ourself to a discrete subset Q C Q. The question is whether 
this restriction will imply a loss of information, that is, whether the set S = {\qk) = l7o-(<jfc))i Qk S 
Q} constitutes a discrete frame itself, with resolution operator 

^=Y1 \lk)iQk\- (8) 

The operator A need not coincide with the original Au. In fact, a continuous tight frame might 
contain discrete non-tight frames, as happens in our case (see later on Sec. H]). 

Let us assume that S generates a discrete frame, that is, there are two positive constants 
< 5 < i? < oo {frame bounds) such that the admissibility condition 

61^1' < E \{Qkm\'<B\m' (9) 

Qk&Q 

holds V'i/' G TC. To discuss the properties of a frame, it is convenient to define the frame (or 
sampling) operator T : H ^ i'^ given by T{ijj) = {{qk\'4'), Qk £ Q}- Then we can write A = T*T, 
and the admissibility condition ([9]) now adopts the form 

bl < T*T < BI, (10) 

where I denotes the identity operator in TC. This implies that A is invertible. If we define the 
dual frame {\q) = A~^\q)}, one can easily prove that the expansion [reconstruction formula) 

= E i^kimk) (11) 

gk&Q 

converges strongly in TC, that is, the expression 

gk&Q gk&Q 

provides a resolution of the identity, where Tf^ = (q'*T)^^T* is the (left) pseudoinverse (see, 
for instance, [17|) of T (see e.g. [15^19] for a proof, where they introduce the dual frame operator 
f = (T+)* instead). 

It is interesting to note that the operator P = TTf^ acting on ("^ is an orthogonal projector 
onto the range of T. 

We shall be mainly interested in cases where there are not enough points to completely 
reconstruct the signal, i.e., under sampling, but a partial reconstruction is still possible. In these 
cases S does not generate a discrete frame, and the resolution operator A would not be invertible. 
But we can construct another operator from T, B = TT* , acting on 

The matrix elements of B are 

Bki = {qk\qi), (13) 

therefore B is the discrete reproducing kernel operator, see eq. ([37|) . If the set S is linearly 
independent, the operator B will be invertible and a (right) pseudoinverse can be constructed 
for T, = T*{TT*)~^ , in such a way that TT^"*" = Ip. As in the previous case there is 
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another operator, Ps = T^T acting on Ti which is an orthogonal projector onto the subspace 
TLs spanned by S. A pseudo-dual frame can be defined as 

\qu) = Y,{^'^)lk\m) (14) 

providing a resolution of the projector Ps, 

T+T = J2 \Qk){qk\ = T*{T+r = ^ \qk){qk\ = Ps (15) 

Using this, a partial reconstruction (an "alias") il^{q) of the signal ip{q) is obtained, 

^(g) = {q\^) = Ekiq)^{qk), (16) 

9fc6Q 

from its samples ^{qk) = {Qk\'^)j through some "sinc-type" kernel (or "Lagrange-like" interpolating 
functions) 

Sfc(^) = (^l^fc) (17) 

fulfilling 'Ek{qi) = 6ki- The alias tjj is the orthogonal projection of tp onto the subspace TCs, that 
is, \ip) = Ps\il^)- The distance from the exact ip to the reconstructed signal is given by the 
error function: 

E4S) = \\i;-i;\\ = y'{tP\I-Psm. (18) 

The two operators A and B are intertwined by the frame operator T, TA = BT. If T is 
invertible, then both A and B are invertible and = Tj^ = T~^. This case corresponds to 
critical sampling, where both operators A and B can be used to fully reconstruct the signal. 

It should be noted that in the case in which there is a finite number of sampling points qt, 
the space £^ should be substituted by C^, and the operator B can be identified with its matrix 
once a basis has been chosen. 



3 Representations of SU{1^1): Discrete Series 

Discrete series representations of SU (1, 1) can be found in the literature (SjH]. Here we shall try 
to summarize what is important for our purposes, in order to keep the article as self-contained 
as possible. 

3.1 Coordinate Systems and Generators 

The group SU{1, 1) consists of all unimodular 2x2 matrices leaving invariant the pseudo- 
Euclidean metric rj = diag(l, —1) and can be parametrized as 

SUih 1) = {UiO = (I I) , Ci, C2 G C : det([/) = |CiP - IC2P = !}• (19) 

The group SU{1, 1) is locally isomorphic to the three-dimensional Lorentz group S0{2, 1) 
(the group of "rotations" of the three-dimensional pseudo-Euclidean space). More precisely 
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SO{2, 1) = SU{1, 1)/Z2, where Z2 = {/, — /} (/ is the 2x2 identity matrix) is the cychc group 
with two elements. The group SU{1, 1) acts on C as 

?7(C) :C^C, z^/ = ^i^^. (20) 

C2Z + Ci 

This action is not transitive, so that C is fohated into three orbits: 

Bi = {z G C : |z| < 1}, C - Bi = G C : \z\ > 1}, Si = {z G C : \z\ = 1}. (21) 

The open unit disk Di may be considered as the stereographical projection of the upper sheet of 
the two-sheet hyperboloid = {{xo,xi,X2) G : Xq — a^f — = 1} onto the complex plane. 
The hyperboloid may be identified with the set of elements of SU{1,1) with Ci = ^0 = 
cosh(r/2) and (2 = +'1x2 = sinh(r/2)e*", r > 0, a G [0,27r[, the stereographical projection 
being given hy z = C2/C1 = tanh(r/2)e*" G Oi. Thus, we could also identify Bi with the coset 
SU{1, 1)/C/(1), where U{1) is the (diagonal) subgroup of the phase e^"^ = Ci/|Ci|- 

The infinitesimal generators of the action (j20p of SU{1^ 1) on C can be written in terms of 
the Lie algebra sn(l, 1) basis elements: 

n CI (X (jj , , 

\AiA/ 

They close the following Lie algebra commutation relations: 

[Ko,K±] = ±K±, [K.,K+]=2Ko. (23) 
It is not difficult to verify that the quadratic (Casimir) operator 

C = K^- {K+K_ + K_K+)/2 (24) 

commutes with every K. 



3.2 Unitary irreducible representations: S'?7(l,l) coherent states 

We are seeking for unitary irreducible representations of SU{1, 1). By Schur's lemma, for any 
irreducible representation of su(l,l), the Casimir operator C must be a multiple of the iden- 
tity /, which we shall set C = s{s — 1)1. Thus, an irreducible representation of SU{1,1) 
is labelled by a single number s (which we shall refer to as the symplectic spin or just "sym- 
pling"). We shall restrict ourselves to discrete series of representations, for which s is half- integer 
s = 1,3/2, 2, 5/2, .... We shall take the orthonormal basis vectors |s, n) in the carrier (Hilbert) 
space Tis to be eigenvectors of Kq: 

Ko\s,n) = {n + s)\s,n). (25) 

From the commutation relations (j23p . we observe that K± act as raising and lowering ladder 
operators, respectively, whose action on the basis vectors proves to be 

K+\s,n) = ^/{n + l){2s + n)\s,n + l), K^\s,n) = ^/n{2s + n - l)\s,n - 1). (26) 

Indeed, it can be easily checked that the action (j25|26|) preserves the commutation relations 
(j23|) : for example: 

[K+,K^]\s,n) = ■■■ = 2{n + s)\s,n) =2Ko\s,n), (27) 
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and so on. From the expression (j26p we deduce that the spectrum of Kq is unbounded from 
above n = 0, 1, 2, . . . , that is, the Hilbert space TCs is infinite-dimensional. 

Any group element f/(C) £ SU{1, 1) can also be written through the exponential map 

U{z,z,^) = e«-f^+-f^-e^'^^«, e = |e|e*^, z = tanh|^|e*^. (28) 

Note that the structure subgroup U{1) C SU{1,1), generated by Kq, stabilizes any basis 
vector up to an overall multiplicative phase factor (a character of ^7(1)), i.e., e^'^^°\s,m) = 
gi(m+s)(/3|^^ Thus, according to the general prescription explained in Sec. [2l letting Q = 
SU (1, 1)/U{1) = Bi and taking the Borel section a : Q ^ G with 17(2;, z) = {z, z, 0), we shall de- 
fine, from now on, families of covariant coherent states mod([/(l), cr) (see [9j). In simple words, 
we shall set (/? = and drop it from the vectors U{z, z, (p)\s, m). 

For any choice of fiducial vector I7) = |s,m) the set of coherent states |z,m) = U{z^z)\'^) 
is overcomplete (for any m) in TLs- We shall use I7) = |s,0) as fiducial vector (i.e., the lowest 
weight vector), so that i^-|7) = and the coherent states 

\z) = U{z,zM = e«-^+-«^-|s,0) =Ms{z,z)e'^+\s,Q), (29) 

are holomorphic (only a function of z), apart from the normalization factor Ng which can be 
determined as follows. Exponentiating the relations (I26p gives 



e 



^■^+|s,0) = |s, 0) + zV2's\s, 1) + ^z2^^2(2s + l)|s, 2) + 

1/2 



= E( r ) '''M^^s{z,z)-'\z). (30) 

n=0 

Then, imposing unitarity, i.e., {z\z) = 1, we arrive at J\fs{z,z) = (1 — l-zp)*. 
The frame 2; G C} is tight in TCs, with resolution of unity 

1= / \z){z\y -2, (31) 

where we denote d^z = (iRe(2;)(iIm(z). Indeed, using (jSOp we have that 



vr Jo^ (1 - zzY vr J^^ y V ™ /V 



(1 - zzf 

2{2s - 1) y + " |s, n)(s, n\ [\l - r^fs-2pn+i^^ 
00 

XI = ^! (32) 



n=0 

where polar coordinates were used at intermediate stage. 
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Using pop , the decomposition of the coherent state \z) over the orthonormal basis {|s,m)} 
gives the irreducible matrix coefficients 

= U'm{z). (33) 

A general sympling s state 

oo 

IV') = ^ am\s,m) (34) 

m=0 

is represented in the present complex characterization by 

oo 

^(z)^(z|^) = J^a^?74(z), (35) 

which is an anti-holomorphic function of The Fourier coefficients a„ can be calculated 
through the following integral formula: 

a„ = {sMi^) = ^ / ^{z^iz) (36) 

Joi (1 - ZZy 

Note that the set of CS {l^;)} is not orthogonal. The CS overlap (or Reproducing Kernel) 
turns out to be 

This quantity will be essential in our sampling procedure on the disk Bi. 

There are other pictures of CS for SU{1, 1) corresponding to other parameterizations like, 
for example, the one that takes Di to the upper complex plane, but we shall not discuss them 
here. 



4 Sampling Theorem and DFT on Di 

Sampling techniques consist in the evaluation of a continuous function ("signal") on a discrete 
set of points and later (fully or partially) recovering the original signal without losing essential 
information in the process, and the criteria to that effect are given by various forms of Sampling 
Theorems. Basically, the density of sampling points must be high enough to ensure the recon- 
struction of the function in arbitrary points with reasonable accuracy. We shall concentrate 
ourselves on sympling s holomorphic functions and sample them at appropriate points. 

In our case, there is a convenient way to select the sampling points in such a way that 
the resolution operator A and/or the reproducing kernel operator B are invertible and explicit 
formulas for their inverses are available. These are given by the set of points uniformly 
distributed on a circumference of radius r: 

Q = {qk = Zk = re^-'^/^, k = 0,1, . . . N - l,r e {0, 1)}, (38) 

^Here we abuse notation when representing the non-analytic functions and simply as Ui^{z) and ^'(z), 
which are indeed (anti-)holomorphic up to the normalizing, non-analytic (real), pre-factor TVs — {l — zz)". Usually, 
this pre- factor is absorbed into the integration measure in pip . 
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which is a discrete subset of the homogeneous space Q = SU{1, 1)/U{1) = Bi, made of the 
roots of , with < r < 1. Denote by 5 = {\zk) , k = 0,1, . . . , N — 1} the subset of coherent 
states associated with the points in Q and by 

7is = Spand^o), \zi), ... , \zn-i)) (39) 

the subspace of 7is spanned by S. For finite we have Tif 7is, so that we cannot reconstruct 
exactly every function ^ £ Tig from N of its samples ^{zk) = (z^lV'); but we shall proof that 
for bandlimited functions 

enf = SpM\s,0),\s,l),...,\s,M)) (40) 
of bandlimit M < oo we can always provide an exact reconstruction formula. 



4.1 Bandlimited Functions 

Theorem 4.1. Given a bandlimited function if) G on the disk Di, of band limit M , with a 
finite expansion 

M 

= ^ am\s,m), (41) 

m=0 

there exists a reconstruction formula of tp 

N-l 

i;{z) = ^k{z)^{zk), (42) 
k=0 

from N > M of its samples ^'(-Zfc) taken at the sampling points in [38\). through a "sinc-type" 
kernel (or "Lagrange-like" interpolating function) given by 

Firstly, we shall introduce some notation and prove some previous lemmas. 

Lemma 4.2. The frame operator T : Ti^^ — > given by T{ip) = {{zk\ip) , Zf^ G Q} [remember 
the construction after Eq. (0)/ is such that the resolution operator A = T*T is diagonal, 
A = diag(Ao, . . . , Am); in the basis (JW ofTL^^, wit^ 

^ N{1 - r2)2- +^ " r2-, m = 0,...,M. (44) 
Hence, A is invertible inTig^ . Therefore, denoting \zk) = A^^\zk), the dual frame, the expression 

N-l N-l 

= ^ \zk){zk\ = ^ \zk){zk\ (45) 

fc=0 k=0 



*The quantities Am are well defined for m € N U {0} and they will be used in the case of band-unlimited 
functions. 
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provides a resolution of the identity in Ti^ . 

Proof. Using (j33p . the matrix elements of T can be written as: 

1 /2 

%n = {Zk\s, n) = i^' - ^) (1 - r'Yr-e-'-^'-/^ = \\I^T^^^ (46) 

where T denotes the Rectangular Fourier matrix (see [1]) given by T^n = -^^e"*^'^'^"/^, A; = 

0, . . . , A/^ — 1, n = 0, . . . , mH Then, the matrix elements of the resolution operator turn out to 
be 

N-\ N-1 
k=0 k=0 

where we have used the well known orthogonality relation for Rectangular Fourier Matrices 
(RFM) (see e.g. Appendix A of [1] for a discussion of some of their properties): 



Af-l 

E 

fc=0 



2ni{n-m)/N\'' ( N , H u = m mod N . 

J =1 0, ifn/mmodA^ ^ = iV5„,^mod7V, (48) 



and the fact that > M. Since A is diagonal with non-zero diagonal elements A^, then it is 
invertible and a dual frame and a (left) pseudoinverse for T can be constructed, Tj^ = A~^T*, 
providing, according to eq. (fT2l) , a resolution of the identity. ■ 

Proof of theorem 14.11 From the resolution of the identity (j45p . any tp G Ti^ can be written 
as 1-0) = ^^~Q ^'(zfc)lzfc), and therefore ^'(z) = (zlV') = J2k=o '^{zk){z\zk). Using that \zk) = 
A~^\zk), we derive that 

M M 

{z\zk) = ^{z\s,m)iA-''UT:, = -7^Y^ X-'/'e'-"^-/''{z\s,m), (49) 

m=0 m=0 

which coincides with Hjt(z) given in eq. (j43p when eq. (jSSp is used. ■ 

Remark 4.3. It is interesting to note that eq. ()42p can be interpreted as a Lagrange-type 
interpolation formula, where the role of Lagrange polynomials are played by the functions 
Lk{z) = Sfc(z), satisfying the "orthogonality relations" Hfc(z/) = P;^, where the operator 
P = TTj^ is an orthogonal projector onto a Af -dimensional subspace of C^, the range of 
T. In the case of critical sampling, A^ = M -|- 1, the usual result Hfc(z;) = 5ik is recovered, but 
for the strict oversampling case, A^ > M + 1, a projector is obtained to account for the fact 
that an arbitrary set of overcomplete data ^{zk)-, k = 0, . . . , N — 1, can be incompatible with 

A reconstruction in terms of the Fourier coefficients can be directly obtained by means of 
the (left) pseudoinverse of the frame operator T: 

Corollary 4.4. The Fourier coefficients am of the expansion lip) = Y^m=QO,m\s,m) of any 
ip G Tig^ can he determined in terms of the data '^{zk) = {zk\'>p) as 



N-1 



am = , m = 0, . . . , M . (50) 



k=0 



^For the sake of briefness, we shall use here the same notation for Rectangular Fourier Matrices as for the 
square ones, namely J^, in the hope that no confusion arises (see Appendix A of [l] for a more precise distinction 
between both cases). 
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Proof. Taking the scalar product with {zk\ in the expression (j4ip of \ip), we arrive at the 
over-determined system of equations 

M 

%mam = ^(^fc), %m = {Zk\s,m), (51) 

m=0 

which can be solved by left multiplying it by the (left) pseudoinverse of T, Tj^ = (^x*T)^^T* = 
A~^T* . Using the expressions of A^^ = diag(AQ ^, A^^, . . . , A^/), given in Lemma HT2l and the 
matrix elements Tkn, given by the formula (j33l) . we arrive at the desired result. ■ 

Remark 4.5. Note that the Fourier coefficients obtained (rectangular) Fourier trans- 

form of the data ^{zk) followed by a rescaling (which can be seen as a filter) by A~^^'^. 



4.2 Band-Unlimited Functions and Undersampling 

When the reconstruction of a band- unlimited function |^) = Yl'n'=o ^nNi ''^) from a finite number 
N of samples is required, we cannot use the results of the previous section since the resolution 
operator A is no longer invertibl^. However, the overlapping kernel operator B turns out to be 
invertible, and a partial reconstruction can be done following the guidelines of the end of Sec. [2] 
(undersampling) . 

Contrary to the case of the sphere [l], where the Hilbert space of functions of spin s, Tig, 
is finite-dimensional, here TCg is infinite-dimensional, and therefore in the partial reconstruction 
of an arbitrary state \^p) = Yl'?^=o'^n\s,n) a considerable error will be committed unless further 
assumptions on the Fourier coefficients a„ are made. Since is normalizable, the Fourier 
coefficients should decrease to zero, thus even if {tp) is not bandlimited, if a„ decrease to zero 
fast enough, it will be "approximately" band limited if the norm of = Yl'^=M+i (^n\s,n) is 
small compared to the norm of {ip), for an appropriately chosen M. Let us formally state these 
ideas. 

Definition 4.6. Let us define by 

M 

Pm = '^\s,m){s,m\ (52) 

m=0 

the projector onto the subspace Ti^ of bandlimited functions of bandlimit M . We shall say that 
a function tp £ Tig on the disk Di, with an infinite expansion 

oo 

IV') = ^a„|s,n), (53) 

n=0 

is a quasi-bandlimited function of bandlimit M and error ej\/ > iff 

(V'|l - Pm\iP) _ E^=M+l . 5 



^While the operator T : Hs has the same expression as in the previous section, the operator A is an 

infinite dimensional matrix given by: Amn = ^]^+pN^]'+qN^n' ^ with m = j + pN and n — j' + qN, that is, is a 
matrix made of N x N diagonal blocks 
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Theorem 4.7. Given a quasi-bandlimited function of bandlim.it M and error cm, V' ^ "^s; there 
exists a partial reconstruction of tp, 

N-l 

i>{z) = ^k{z)-^{zk), (55) 

fc=0 

from N of its samples ^(zk), taken at the sampling points in i38\). up to an error il8\) 

^ < 4. ^ ^; - ') "> + o(.-) ,5.) 

provided that N > M . The sinc-type kernel |j?D now adopts the following form: 

^^<^' = K& E^7'EW(«.T-F-", (57) 

w/iere 

oo 

= J^Aj+gjv, j = 0,...,iV-l, (58) 

are the eigenvalues of the discrete reproducing kernel operator B = TT* [defined in with 
matrix elements Bki = {zk\zi)] and Xn is given by [44^ , but now for n = 0, 1, 2, . . . . 

Note that the quadratic error ([56|) approaches as r — > and/or N ^ oo. Before tackhng 
the proof of this theorem, we shall introduce some notation and prove some previous lemmas. 

Lemma 4.8. The pseudo-frame operator T : Tig given by T{ijj) = {{zk\tlj),Zk G Q} 

[remember the construction after Eq. (0)/ is such that the overlapping kernel operator B = TT* 
is an N X N Hermitian positive definite invertible matrix, admitting the eig en- decomposition 
B = TT)T* , where D = diag(Ao, • • • , Aat-i) is a diagonal matrix with Xj given by ^58\) and T is 
the standard Fourier matrix. 

Proof. Let us see that B is diagonalizable and its eigenvalues are given by the expression 
()58p . which actually shows that all are strictly positive and hence B is invertible. This can be 
done by taking advantage of the circulant structure of B (see e.g. Appendix B in [1]). Actually, 
using the expression of the CS overlap (f37|) . we have: 

/ l-r2 

Bki = {zk\zi) = [j^^j^-(pj7^ ) ^ Ci_k, (59) 

where the circulant structure becomes apparent. The eigenvalues of B are easily computed by 
the formula: 

N-l 
n,m=0 

Expanding 

Ci = {l- r^f' f2 ^ J " ^) r2'?e2-ViV^ (61) 

and using the general orthogonality relation for Rectangular Fourier Matrices (j48|) . we arrive at 
(j58|) . It is evident that A^ > 0, V/c = 0, 1, . . . , — 1, so that B is invertible. ■ 
Following the general guidelines of Sec. 2, we now introduce the projector P^: 
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Lemma 4.9. Under the conditions of the previous Lemma, the set{\zk) = Z]/=o {^~^)ik\zi) , k = 
0, . . . , — 1} constitutes a dual pseudo-frame for S, the operator Pg = T^T is an orthogonal 
projector onto the subspace TLf , where = T*B^^ is a (right) pseudoinverse for T , and 

N-l N-1 

\zk){zk\ = ^ \zk){zk\ = Ps (62) 
fc=0 fc=o 

provides a resolution of the projector Ps, whose matrix elements in the orthonormal base 125\) 
of Tig exhibit a structure of diagonal N x N blocks: 

{s,m\Ps\s,n) = Pmn{r,N) = {\mK)^^'^Ki]nodN^n,mraodN-, m,n = 0,... ,M, (63) 

with An and A„ given by [58\) . 

Proof. If we define T^^ = T*B~^ it is easy to check tliat TT^^ = In is tlie identity in C^. In 
the same way, Ps = T^'^T is a projector since Pj = T^TT^T = T^T = Ps and it is orthogonal 
Pg = {T*B~^T)* = T*B~^T = Ps since B is self-adjoint. The resolution of the projector is 
provided by Eq. (jlSp . Its matrix elements can be calculated through: 

N-l 

Pmn{r,N) = T.L{B-^)ik'Tkn- (64) 

k,l=Q 

The inverse of B can be obtained through the eigen-decomposition: 

N-l 

{B-')ik = {Tb-^r')ik = ^ E A-ie2-^(^-')/^. (65) 

i=o 

Combining this expression with (j46p . and using the general orthogonality relation (j48p for RFM, 
we finally arrive at (|63p. ■ 

The matrix elements (|63p will be useful when computing the error function (jl8l56p for a 
general quasi-bandlimited function (I53ti54p . We are interested in their asymptotic behavior for 
large (large number of samples) and small r. 

Lemma 4.10. Denoting m = j +pN and n = j + qN , with j = 0, . . . , N — 1 and p,q = 0,1, ... , 
we have the following asymptotic behaviour for the matrix elements (EJ) of the projector Ps: 

P^„(r,iV)=0(r(^'+«)^). (66) 



Proof. Let's define: 

u=l \ n j \ n ) 

Using ([63]) we have: 



^2 /2s+i+piV-lNl/2/2s+i+(?Af-l\l/2 
p AT\ _ [^j+pN^j+qN) _ \ j+pN J \ j+gN J Jp+q)N 



E^=oA,+.iv (2^7-1) l + e,{r,N) 

(2s+j+pN-l\ 1/2 /2s+j+qN-l\ 1/2 
I j+pN 2 1 

m 

which gives the announced asymptotic behaviour (j66l) for small r and/or large A^.l 



J+P^ ^^J+g^ ^ ^(p+a)jV(X - 0(r2^)), 
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Lemma 4.11. The functions |g7| j are decreasing sequences of n, that is: 
en{r, N) < em{r, N) <^ n > m, n,m = 0, . . . , N — 1. 



Proof. The sequence Enir^N^ is decreasing in n since the quotient of binomial coefficients 
^^~n+uN^^) I (^^"n'''") (|671) is decreasing in n for any u G N, as can be seen by direct compu- 
tation. ■ 

Now we are in conditions to prove our main theorem in this section. 

Proof of Theorem (I4.7p 

Introducing the expression ()65p in (jl4p (with = z^) and noting that 

oo 

{z\zi) = ^^{z\s,n){s,n\zi), (70) 

n=0 

we arrive at ([57|) after using p!7|) . ([33]) and the orthogonaUty relation (08]). 

Now it remains to prove the asymptotic behaviour (|56p for the error. Decomposing \'\\)) in 
terms of |V'm) = -Pa/IV') and IV'm) = ~ -Pm)|V')) we can write 

4(r,iV) = (V.jv/|V'M) + (V'i^|V'M)-(V'Af|^'5|V'M)-(V'Ml^'5|V'M) 

-2i?e(V^M|P5|V'M) (71) 

For simplicity, let us assume that iV = M + 1. Using that ■0 is a quasi-band limited function 
and the results of Lemma 14.101 and 14.11] we have that 

4(r,iV) < el,m) + (1 - m) - k^iAPsM^ii) - IRc^i^um^ii) (72) 

The third term to the r.h.s. is positive and small (of 0(r^^), according to Lemma l4.10p . 
and can therefore be neglected. The last term has no definite sign, and can be seen to be of 
0{r^), so it must be taken into account. Bounding it in absolute value and using repeatedly 
the Cauchy-Schwarz inequality we have that: 



Using (|67p we arrive to the desired result. ■ 
Remark 4.12. Given e > cm, an upper bound for r and a lower bound for can be given in 

order to have |^^|^^ < e^, by solving numerically ([75]) < e^. Equation ([^5]) can be used to 
obtain an analytic estimate of an upper bound for r, which turns to be: 



i_ 

N 



(74) 
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Corollary 4.13. (Discrete Fourier Transform) The Fourier coefficients an of the expansion 
can be approximated by the discrete Fourier transform on the hyperboloid: 

.1/2 , N-l 
^n 



1 J2 e2™'=/^xl/(z,), (75) 



Proof. 

an = {s,n\7p) = {s,n\Ps\ij) = ^{s,n\zk)'i'{zk) = ^ {s,n\zi){B-^)ki'i>{zk) (76) 

fc=0 k,l=0 

and this provides the desired result once (06]), ([65|) and the orthogonahty relation (|18]) are used. 



N-l N~l 



Remark 4.14. The discrete Fourier transform on the hyperboloid resembles the traditional dis- 

1/2 ^ 

Crete Fourier transform up to a normalizing factor An /AnmodiV) which can be seen as a filter. 
Remark 4.15. The definition of the Fourier coefficients entails a kind of "periodization" of 
the original a„ in the sense that 



q=0 

which implies 



On = ^ A^^Aj- ^Aj^^^^^Oj+gAT, j = nmodN, (77) 

q=0 



,_l/2~ \-l/2 - ^ - l^n+pN . 

fln — \+pN^n+pN ^ On+piV — y — U"j 

We could think that, for the case eA/ = 0, we should recover the results of Section BTl but 
we shall see that this is not the case. Before, a process of truncation and rescaling of {tp) in 
(j55p is necessary to recover the reconstruction formula ()42p for strict bandlimited functions ()4ip . 
Indeed, the truncation operation 

M 

\i^M) = Pm\iP) = ^ an\s,n) (79) 

m=0 

followed by a rescaling of the Fourier coefficients 

M , 

IV'm) = R\^m) = V ^an|s, n) (80) 

n '^'l 

renders the reconstruction formula for i^^j{z) = (zli^P/v/lV') to the expression (j42p . 

This fact suggests us an alternative approach to the sampling of quasi-bandlimited functions 
■0 for small em, which will turn out more convenient in a certain limit. Actually, for cm << 1 
we have 

\\i^-PMn' = eii\m'«\m', m 

so that, the reconstruction formula for ipM = Pm'4' would give a good approximation of 
■i/;. The problem is that, in general, the data ^(zf^) = (-ZfciV') for ip and the data ^M^Zk) = 
(zfelP/Vf IV') for tpM are different unless {zk\PM = VA; = 0, . . . , A'' — 1, which is equivalent to 
{zk\PM\zk) = l,yk = 0, . . . , N — 1. The following proposition studies the conditions under which 
such requirement is satisfied. 
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Proposition 4.16. For large sympling s ^ oo and large band limit M oo, the diagonal 
matrix elements of Pm in Tlf have the following asymptotic behaviour: 

PMr) - (..iPMk.) - { ; i; IX^Jl (82) 



where 



(i + \r^)-'/' (83) 



M 

denotes a critical radius. For M >> 2s >> 1 we have < 1. 
Proof. Using the expression (|46|) we have 

M /9 _i_ _1\ 

P!,{r) = {z,\PM\zk) = E %A = (1 - E ( r 

171=0 m=0 

Denoting by p = we can compute 



(85) 



We identify here the Binomial distribution B{2s — 1 + M,p) (up to a factor 2s + M), which has 
a maximum (as a function of p) for pc = 1/(1 + ^fy^)- Using the Central Limit Theorem for 
2s — 1 + M ^ oo and the representation of the Dirac delta function as the limit of a normal 
distribution, we identify (j82l) as a Heaviside-type function, concluding the desired resultH 



5 Conclusions and Outlook 



We have proved sampling theorems and provided DFT for holomorphic functions on Bi carry- 
ing a unitary irreducible representation of SU{1, 1) of sympling s. To accomplish our objective, 
we used the machinery of Coherent States and discrete frames, and benefit from the theory of 
Circulant Matrices and Rectangular Fourier Matrices to explicitly invert resolution and repro- 
ducing kernel operators. We also paved the way for more general coset spaces Q = G/H and 
their discretizations. 

Heisenberg-Weyl (and Newton-Hooke) CS could be seen as a zero curvature limit k ^ (and 
large s) of SU{2) (positive curvature) and SU{1, 1) (negative curvature) CS, a unified treatment 
of sampling for the three cases being possible. This is left for future work [18j. 
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